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^ : Abstract 
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It is pointed out that the recent measurement of the angle 7 of the unitarity triangle, 

_^ , providing irrefutable evidence for a complex Cabibbo-Kobayashi-Maskawa (CKM) 

f— ^ ' matrix, presents a great challenge for supersymmetric models with spontaneous CP 

\^ • violation. We construct a new minimal extension of the minimal supersymmetric 

standard model (MSSM), with spontaneous CP breaking, which leads to a com- 
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T^ I plex CKM matrix, thus conforming to present experimental data. This is achieved 

II ■ through the introduction of two singlet chiral superfields and a vector-like quark chi- 

Oh. ral superfield which mixes with the standard quarks. A Z3 symmetry is introduced 

(— I I in order to have a potential solution to the strong CP problem. 
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1 Introduction 

Four decades after the experimental discovery of CP violation, the origin of CP 
breaking remains one of the fundamental open questions in particle physics. 
There are two basic scenarios for CP violation in the framework of gauge the- 
ories, namely explicit CP breaking at the Lagrangian level and spontaneous 
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CP breaking by the vacuum (SCPV) [1,2]. In the Standard Model (SM), one 
assumes that CP is exphcitly broken at the Lagrangian level through the intro- 
duction of complex Yukawa couplings leading to CP violation in the charged- 
weak interactions, parametrised by a complex Cabibbo-Kobayashi-Maskawa 
(CKM) matrix [3]. 

In supersymmetric (SUSY) extensions of the SM there are additional sources 
of explicit CP violation, arising from complex soft SUSY breaking terms as 
well as from the complex SUSY conserving yU parameter. The SUSY phases 
also generate large contributions to the electric dipole moments (EDMs) of 
the electron, neutron and mercury atom. The non-observation of the EDMs 
imposes strong constraints on the SUSY phases, forcing them to be very small. 
This is the so-called SUSY CP problem and many solutions have been pro- 
posed to overcome it [4,5]. 

SCPV is an attractive approach to the SUSY CP problem, since all the cou- 
plings of the Lagrangian are real, due to the imposition of CP invariance at the 
Lagrangian level. The only source of CP violation are the vacuum phases [1,2]. 
SCPV also provides an appealing solution to the strong CP problem, since in 
this case one may have naturally vanishing O at tree level and calculable at 
higher orders [6,7,8,9,10,11]. 

In this paper, we construct a minimal supersymmetric extension of the SM, 
with spontaneous CP breaking, which is in agreement with all the present 
data, provided by both BaBar and Belle collaborations [12,13]. The question of 
compatibility of spontaneous CP breaking in supersymmetric (SUSY) models 
with recent CP data is highly non-trivial, for the following reason. One of 
the salient features of most SUSY models with spontaneous CP breaking is 
the fact that they lead to a real CKM matrix, since the phases in quark 
mass matrices can be removed by rephasing of right-handed quark fields, in 
a manner entirely analogous to the one encountered in models with natural 
flavour conservation in the Higgs sector [2]. Until recently, SUSY models with 
spontaneous CP breaking and a real CKM matrix [14,15] could be viable since 
both Sk and aj/^^s could be generated by supersymmetric contributions to 
j{0 _ j{0 g^]^(j j^^ _ j^^ mixing, respectively [16]. The novel experimental input 
is the recent measurement of the angle 7 of the unitarity triangle [12,13]. A 
recent analysis [17] of the present experimental data provides clear evidence 
for a complex CKM matrix even allowing for the presence of new physics 
contributions to Ek, clj/iPKsi ^^b^ and AMb^- These experimental findings 
lead to the question of whether it is possible to have a SUSY extension of 
the SM with spontaneous CP violation and a complex CKM matrix. In this 
paper we will show that the answer to the above question is in the affirmative. 
Indeed we construct a minimal supersymmetric extension of the SM (MSSM) 
with spontaneous CP breaking which leads to a complex CKM matrix, thus 
avoiding confiict with recent data. The crucial point is the introduction of 
at least one Q = —1/3 vector-like quark which mixes with standard quarks 
and leads to a non-trivial phase in the 3x3 effective CKM matrix connecting 



standard quarks. The existence of such matter states naturally arises in various 
extensions of the SM, such as in Eq grand unified theories as well as in models 
with extra-dimensions, including some of the superstring inspired models. 
The model presented here has another advantage with respect to previously 
suggested SUSY models with SCPV. This stems from the fact that physical 
phases only arise in the vacuum expectation values (VEVs) of Higgs singlet 
phases. As these singlets only couple to the quarks and squarks, they will not 
show up in the neutralino and chargino sectors, making the solution of the 
SUSY CP problem easier in this context. 

2 The Model 

We consider a simple extension of the MSSM where in addition to the usual 
superfield content, we introduce an isosinglet vector-like down quark, D4, D^, 
and in the Higgs sector, two singlet chiral superfields Si, S2 which have, under 
SU{3)c X SU{2)i^ X U{1)y, the following quantum numbers 

D, : (3, l)_,/3 , D^^ : (3, 1),/, , 81,82: (1, l)o • (1) 

We will show that the extended Higgs structure can spontaneously break CP 
through phases in the VEVs of the singlet scalar components Si, 82- Further- 
more, these phases induce a non-vanishing phase in the CKM matrix, which 
is not suppressed by the small ratios \{hu,d)/{si,2)\ ^ 1- We emphasise that 
small values of the later ratios are crucial in order to naturally suppress flavour- 
changing neutral currents (FCNC), but also the one-loop flnite contributions 
to the parameter 6 associated with strong CP violation [18]. 
Since we want to achieve spontaneous CP breaking, we impose CP invariance 
at the Lagrangian level, which implies that all the parameters appearing in 
the superpotential and in the soft breaking terms are real. Moreover, we intro- 
duce a discrete Z3 symmetry, under which the standard superflelds transform 
trivially, while the new chiral superflelds have the following Z^ assignment 

{D^',8,,82} ^e'-^^ D4-e^\ (2) 

The role of the Z3 symmetry is to forbid quark bare mass terms of the form 
dihd^R, thus leading to the vanishing of at the tree level. The most general 
renormalisable gauge invariant superpotential which is compatible with R- 
parity invariance and the Z3 symmetry reads as 

W=e,,{ Y:^ Q1 Hi Uf + F/ Q', H^, D^ + Y;^ L\ H^, Ef + /x H^ H^, ) 

+ f 8i Di D^ + g' 82 Di D^ + Md D^ D^ (3) 

+ All ^l + A12 Oi 02 + A2I Ol 1S2 + A22 '-'2 5 

where the indices a,h = 1,2 are 8U{2) indices and e is a antisymmetric 2x2 



matrix, with £12 = 1. We have assumed that the singlets superfields Si, S2 
are even and the isosinglet superfields D4, D^ are odd under the extended R- 
parity. The matrices F^ , F^ , Y^., the vectors / , g and Mj:, generate the masses 
of quarks (including the vector-like one) and leptons. 

In addition to the superpotential given by Eq.(3), we have to specify the 
explicit soft-breaking terms, which read as 

+ (M|).,efe^- + MlJJ: + [Ml^hrK+Ml^hl*hi) 

+ Ml si si + Ml si s, + e,, (A^ Y^^ q^ u'r hi + A'i Y^ q\ d^ K + 

+ AiY:^~l'^e^ K + B ^^Kh', + h.c. 

+ {Afi /* si di d^ + Agi g^ S2 d^ d^ + h.c. 

+ (Aau All si + Axr, A12 si S2 + Aa2i A21 Si s^ + Ax^^ A22 ^2 + h.c. 



(4) 



- (Mi BB + M2 WW + M3GG + h.c. 



In the above equation, we take the soft-SUSY breaking coefficients Mq , B, 
Au, ... and Mi^2,3 as free parameters at the weak scale, and choose them real in 
order to respect CP invariance of the Lagrangian. 



3 Analysis and results 

3.1 Minimisation of the neutral scalar potential 

After spontaneous electroweak symmetry breaking, in the neutral Higgs sector, 
we have verified that there is a region of parameter space where minimum of 
the Higgs potential is at: 

(M=r I , {hu)=\ ° , {si) = Vie"f\ {s,)=V2e''^\ (5) 

where Vd, v^, Vi are real, positive and 6, 0i and 02 are physically meaningful 
phases. It will be shown that the minimum of the potential is at ^ = 0. On the 
contrary, the phases 0i, 02 are in general non- vanishing and lead to SCPV. 
Furthermore, we will show that these phases lead to an effective down quark 
mass matrix capable of generating a CP violating phase in the 3x3 sector of 
the CKM matrix. 



From the superpotential and soft SUSY breaking terms, Eq.(3) and Eq.(4), 
we derive the following CP-invariant neutral scalar potential: 

Ke... = K= + Ml V^ + Ml V^ 

9 A^^ + 1 \l,) Vt + 4 (\l, + A^^) V^ Vl + (9 \l^ + 1 A^,) V^ 

+ 4 Vi \/2 (3 An Ai2 Vl + A21 (3 A22 V^ + A12 {Vl + V^))) cos (0i - ^2) ,^, 

(6) 

+ 6 All A21 Vl Vl cos (2 (01 - 02)) + 6 A12 A22 Vl V^ cos (2 (0i - 02)) 
+ 2 A),,, A12 V;' V2 COS (2 01 + 02) + 2 Aa,, A21 Vi V^ cos (0i + 2 02) 
+ 2 Aa,, All V^ cos (3 0i) + 2 A,, A22 V2' cos (3 02) . 
where the MSSM part is 

K^r? = (Mi, + /.^) .^ + (Mi„ + ^^) .^ + i (/ + ,' 

— 2 _B/i VuVd cos ^ . 



"^^-^--^)' 



(7) 



From Eq. (7) it is clear that the minimisation equations will require ^ = 0, 
for /5/i positive. Note that /9/i positive is required, as in the MSSM, in order 
to obtain positive squared masses for the physical CP-odd state contained 
in hui hd- This is to be expected as the singlets do not mix with the MSSM 
doublets and it is well known that the MSSM does not to lead spontaneous 
CP violation [2,19]. Setting ^ = we still have six minimisation equations for 
the four VEVs and two phases. The minimisation of this potential is quite 
involved. We have done it numerically, following the procedure of Ref . [20] . To 
explain it briefly, we write the minimisation equations in the form: 

t^*' neutral c^ ( Ti f2 , 2\ U ^ U ( 2 2\ r, D 

g^^ = 2 [Mh^ + f^)Vd [v^ - V^j Vd-2 B^Vu , 

ni/ „2 I „/2 

l^*' neutral r) / » ^2 |2\ , iJ ^ iJ f 2 2\ r, D 

—^ — = 2 [Mh^ +f^ )v^ + [v^ - v^j v^-2 BfiVd , 

8V 

-^ = 2 MlV^ + Fi(A,,, A,^^,V„ 0,) , (8) 

8V 

— ' 2MlV, + F,{X,,,A,^^,V.,<j>,), 



dV, 



2 



^^ = 6AA,,AiiV;=^sin0i + F3(A,„AA,,,l^„0.), 

0(pi 

dV 
aT"' = 6 Aa22 A22V'2^ sin 02 + Fi{\ij, Ax , 14, 0i) . 

1702 

These equations can be analytically solved for the quadratic soft masses [20] 



and for the Aa^ and Ax^^ parameters [15]. The procedure is therefore to scan 
over all the VEVs, phases and remaining parameters, and obtain Mfj , Mfj^, 
Mg_^, Mg^, A^ii, Ax22 from Eq. (8). We are then certain that the extrema 
equations are satisfied. To assure that we are at a minimum we compute 
the eigenvalues of the neutral Higgs boson mass matrices and require that 
they are all positive, except for the Goldstone bosons that, of course, should 
remain massless. If this is verified we are at a minimum. Finally we check if this 
minimum is deeper than the minimum that does not violate CP. Following this 
procedure we easily obtained a very large number of good solutions, which are 
true minima of the potential and do violate CP. Next we analyse the question 
of CP violation in the quark sector. 

3.2 The quark mass matrices 

Upon electroweak gauge symmetry breaking, quark mass matrices are gener- 
ated through the Yukawa terms, Yu ,Yd, f ,g, and the mass term Md. Note 
that Yu , Yd and Md are real. The couplings / and g are also real, but the mass 
terms generated by them are in general complex due to the phases (pi = arg(sj). 
Thus, the quark mass matrix are encoded in the following Lagrangian: 

-^MAss = uJ{m^)ijU^ + d^{Md)apd^ + h.c. , (9) 

where i, j = 1, 2, 3, a = 1, 2, 3, 4, and the matrices m„ and M^ are given 

(nid 



rriu = VuYu, Md 



(10) 
\M Md, 



with rud = VdYd, and Mj = fiViC^'^^ + giV2e^'^'^. In order to simplify our 
analysis, and without loss of generality, we perform an orthogonal weak basis 
transformation that leaves the matrix Yu in a diagonal form. The complex 
elements Mi can be written in terms of moduli and phases as 



M = [\Mi\ e"^\ IM2I e*'^% IM3I e*'^^ j ^ (n) 

where these moduli and phases are given by 

l^il' = f'M + alV:^ + V^9^VlV2 COS(01 - ^2) , (12) 

and 

/iVi sin 01 + 5(^1/2 sin 02 .,„^ 

/jl/lCOS0i +5'iV/2COS02 

respectively. Moreover, the phase (pi can be set to zero, without loss of gen- 
erality, by re-definition of the vector-like quark fields. The phases (^2 and 793 
are the only source of CP violating that appear encoded on the 3x4 CKM 
matrix. In what follows, we assume rrid to be a 3 x 3 real symmetric matrix. 



Since the mass terms M , Ma are A J = 0, they can be much larger than Vu 
and Vd- 

As we have mentioned in the introduction, the strong CP problem is auto- 
matically solved at tree level. The parameter G associated with strong CP 
problem can be written as 9 = 6qcd + ©qfd, where 6qcd is the coefficient of 
g'^FF/32TT'^ which vanishes since CP symmetry is imposed at the Lagrangian 
level. The parameter 6qfo, at tree level, is given by 

Gqpo = arg [det(m„) det(Md)] , (14) 

since the matrix rriu and the determinant of the matrix M^ from Eq.(lO) are 
real, 

det (Md) = Md det (m^) , (15) 

it follows that 

e,... = 0, (16) 

thus providing a strong CP problem solution of the Barr- Nelson type [10,11,21]. 
Note that the determinant of Eq.(15) is real due to the Z^ symmetry, which 
forbids terms like Qi H^D^. 

In order to determine the left-handed unitary rotations which encodes the 
quark-isosinglet mixings and the CP violation information one has to diago- 
nalise the hermitian matrix H = Md M\ (4 x 4) written as 

^ = I _, \. (17) 

Mrud M 

where M = |Mip + IM2P + jMsp + M|,. The matrix H can be diagonahsed 
by the following unitary matrix [18] 



U 



leading to: 

U^ HU = \ ^ , (19) 

VO Ml] 

where V^ is diagonal formed with the light down quark squared masses. From 




Eq.(18) and Eq.(19) one derives: 



J- rudM'fMml f^ H \ ^ 



An effective 3x3 light down quark mass matrix can be easily obtained in the 
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limit where K is an unitary matrix and H <^ M [18,22]. One obtains: 

meflfm^jj ^rud I 1 - ^ I "^^ ' (21) 

while the mass of the down-type quark isosinglet is given in an excellent ap- 
proximation by 

Md^M. (22) 

Note that the matrix K is almost unitary since from Eq.(19) and exact uni- 
tarity of U, it follows that 

K^K = 1-S^S, (23) 

with S ~ 0{m(i / M) [18]. Although, the matrix irid is real, the phases in 
M'I'M in Eq.(21) are sufficient to generate a complex K . The full CKM matrix, 
which now has dimension 3x4, comes directly from the matrix U and is given 
by: 

{Vc^mU =(k r) , (24) 



ta 



which contributes to the electroweak Lagrangian as: 

^W,Z = —y= {UL^Y {VcKM)i^dLaW+ + h.C.) + 
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2 cos ^w 
where 



ul a^UL i - di al'^X^p dip -2 sin 9y/ J^^ 1 Z^ , 



(25) 




Xa^^{vL:VcKu)^=S^^- , (26) 



'a/3 

and the electromagnetic current, Jem 5 is given by 

2 1 - 

Jem = i^ua^Ui - -d^-f^d^ . (27) 

The 3x4 CKM matrix from Eq.(24) encodes three independent phases which 
violate CP. However in the limit where one neglects the small mixings between 
the physical heavy quark ^4 and the light quarks {i.e. the limit where K 
becomes unitary) only one phase appears in K [23]. 

3.3 Numerical Example 

In this section we present one concrete numerical example which leads to 
values of quark masses and mixings in agreement with experiments. It is clear 
from the effective mass matrix presented in Eq.(21), that the light quarks and 
their mixing do not depend on the overall scale, which only rescales the heavy 
vector-like down quark mass. We have verified numerically that this feature 



still holds, at a very high precision, provided Mo > 1 TeV, which enforces the 
goodness of the effective matrix given in Eq.(21). We have then parametrised 
the input values Mj as a function of an overall scale M^ = 1 TeV. Moreover, 
we fixed the down quark masses to the following values: 



rud = 4.38 MeV , m, = 94.6 MeV , mt = 3. 1 1 GeV . 



f28) 



We have verified, using Eq.(20) that there are solutions to the matrix rrid 
leading to the effective 3x3 CKM matrix K consistent with experimental 
data. This was done by varying the ratios Mi/Mo and the others parameters 
in a consistent way with Eq.(20). A particular example is obtained with the 
following values for the ratios Mi/Mo'- 



Ml 



0.5, 



Mo 



0.8, 



M, 



3.8, 



Md Md Md 

and the following values for the phases ipi (in radians): 

(P2 = 1-1, V^3 = 1-7. 

We considered the moduli of matrix K elements as 

^ 0.97 0.22 0.0038^ 
= \K\ = 0.22 0.97 0.04 

y0.0086 0.04 0.9992/ 



(29) 



(30) 



K 



3x3 



(31) 



that are within the range allowed by experiments [24]. The relevant CP vio- 
lating weak basis invariants that encodes phases content in K are 



J = Xm{Ki2Kl^Kl^K2z) = 3.14 x 10" 



and 



(3 = arg(-ir2ii^2*3^3i^33) = 23.8° 
which corresponds to: 

sin 2/3 = 0.739, 



(32) 

(33) 
(34) 



in agreement with the most recent experimental data [24]. Note that the quan- 
tities J and sin 2(3 are only relevant to describe CP violation in the quark sector 
when K is nearly unitarity. The deviation from unitarity of the matrix K is 
given by 

1 - K^K = S^S c:^ O (lO-^) , (35) 



which confirms our previous assertions about the almost unitarity of K. Thus, 
the values of the elements of m^ obtained are: 



rrid 



/O.OI 0.025 0.02\ 

0.025 0.13 0.47 

V0.02 0.47 9.67/ 



(36) 



Once the matrix irid is determined, the mass of the heavy vector-like down 
quark state can be calculated by using Eqs. (17-19), giving the following result: 

M^~M = 4.0TeV, (37) 

which enforces the excellent approximation given in Eq.(22). 
Finally we have to verify if the inputs are compatible with the minimisation 
of the potential. The only parameters that enter the quark sector and that 
are relevant for the minimisation of the potential are the VEVs Vi, V2 and 
the phases 0i and 02- As can be seen from Eq. (12) and Eq. (13) there is a 
significant freedom in obtaining the values of Vi and 0j from Mj and (fi. We 
have numerically verified that there is a large number of possible choices of 
parameters that give a successful minimisation for a given choice of Mj and 
(/9j, thus implying that our solutions are fully consistent. 

4 Discussion and Conclusions 

We have presented a new minimal supersymmetric extension of the SM, where 
CP is spontaneously broken. The novel feature of the model is the fact the 
CKM matrix is non-trivially complex, in contrast with previous SUSY models 
with SCPV. 

It was emphasised that having a complex CKM matrix is crucial, in view of 
the recent measurement of the angle 7 of the unitarity triangle. Prior to this 
important experimental result, it was possible to have viable SUSY extensions 
of SM, with SCPV and a real CKM matrix, where CP violation in the kaon and 
B-sectors was entirely generated by new SUSY contributions to K^ — K^ and 
B — B mixings. This class of models are no longer viable. In order to generate 
a non-trivial CP-violating phase in the CKM matrix, we have introduced a 
Q = —1/3 heavy isosinglet quark. The mass terms mixing this heavy quark 
with the standard quarks are responsible for the generation of a complex 
effective 3x3 CKM matrix, which in turn leads to a complex CKM matrix. 
Due to the appearance of naturally suppressed FCNC at tree level, in the 
down-quark sector, there are in the model new contributions [25] to K^ — K^, 
Bd—Bd and B^ — Bg mixings. At present, all experimental data on CP violation 
as well as on rare K and B decays is in impressive agreement with the SM. It 
seems by now clear that the SM and its CKM mechanism of mixing and CP 



10 



violation give the dominant contributions to the physical quantities entering 
the standard unitarity triangle. However, there is still room for New Physics, 
that can even give a dominant contribution to physical quantities which do 
not enter directly in the standard unitarity triangle, an example being the 
invariant phase x — ^T^&{~'VtsVflV*gVcb) ■ This phase has not been measured 
yet, but it will be measured at LHCb. 

We have also emphasised that the class of models presented here, have the 
potential of alleviating the so-called SUSY CP-problem. In particular, we have 
pointed out that since CP violation arises entirely from VEVs of the Higgs 
singlets, they will not show up in the neutralino and chargino sectors. 
A complete analysis of the phenomenological implications of this model, in- 
cluding the evaluation of EDMs is outside the scope of this paper and it will 
be presented elsewhere. 

The fact that a minimal realistic extension of the MSSM, capable of generat- 
ing spontaneous CP violation, requires the introduction of vector-like quarks, 
provides further motivation for the search of these heavy fermions which arise 
in a variety of extensions of the SM. 
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